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Abstract 

Polarization  characteristics  of  metal- supported  solid  oxide  fuel  cells,  fabricated  by  pulsed  laser  deposition  and  suspension  plasma  spray,  were 
investigated.  With  the  electrochemical  impedance  spectroscopy  and  current-voltage  polarization  curves,  the  metal- supported  SOFCs  were  compared 
and  analyzed  in  terms  of  their  cathode  exchange  current  density,  polarization  loss  and  maximum  power  density  over  the  temperature  range  of 
400-600  °C.  The  electrochemical  mechanism  for  the  linear  polarization  characteristics  observed  from  the  experimental  data  was  addressed.  Results 
from  experiments  and  simulation  indicate  that  fabrication  processes  and  operation  temperatures  play  an  important  role  in  the  electrochemical 
mechanism  for  linear  polarization  characteristics  of  metal- supported  SOFCs. 
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1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  expected  to  be  highly 
favoured  for  stationary  power  generation  due  to  their  high  effi¬ 
ciency,  low/zero  pollution  and  fuel  flexibility.  However,  high 
cost  and  low  durability/reliability  are  major  barriers  in  the  path 
towards  SOFC  commercialization.  It  has  been  recognized  that 
both  these  barriers  are  primarily  due  to  the  high  temperatures  of 
operation  [1]. 

In  an  effort  to  overcome  these  barriers,  recently,  SOFCs 
with  reduced  operation  temperature  (500-850  °C)  have  received 
increasing  interest.  SOFCs  operating  at  lower  temperatures 
exhibit  numerous  advantages  such  as  a  wider  choice  of  low-cost 
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component  materials,  improved  stability,  and  increased  flexi¬ 
bility  in  the  structure  design  [2].  As  a  representative  of  SOFC 
family,  metal- supported  SOFCs  are  expected  to  be  more  com¬ 
petitive  than  traditional  energy  conversion  technologies  due  to 
the  low  cost,  high  strength,  better  workability,  good  thermal  con¬ 
ductivity,  and  quicker  start-up  [3].  However,  their  performance 
is  still  relatively  low  at  this  current  technology  stage. 

With  respect  to  the  development  of  reduced  temperature 
SOFC  technology,  reducing  electrolyte  thickness  and  develop¬ 
ing  alternative  materials  with  high-ion  conductivity  at  reduced 
temperature  are  two  major  approaches.  So  far,  various  process¬ 
ing  techniques,  including  spray  pyrolysis  (SP)  [4,5],  pulsed  laser 
deposition  (PLD)  [6-8],  and  suspension  plasma  spray  (SPS) 
[9,10],  as  well  as  spray  and  spin  coating  [11-13],  have  been 
extensively  explored.  These  processing  techniques  have  a  sig¬ 
nificant  effect  on  the  SOFC  interfacial  morphology  that  in  turn 
affects  or  determines  the  SOFC  performance. 

For  performance  improvement  and  component  optimiza¬ 
tion  of  SOFCs,  fundamental  understanding  of  the  correlation 
between  fabrication  processes  and  fuel  cell  reaction  kinetics 
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including  analysis  of  polarization  characteristics  seems  to  be 
necessary  [14-19].  However,  for  metal- supported  SOFCs,  there 
is  very  little  literature  focusing  on  the  fundamental  under¬ 
standing,  in  particular,  through  polarization  analysis  assisted  by 
electrochemical  impedance  spectroscopy  (EIS)  [20]. 

In  our  previous  work  [8],  EIS  and  electrochemical  polariza¬ 
tion  were  employed  to  investigate  the  performance  and  aging 
characteristics  for  metal- supported  SOFCs,  named  Cell  #1,  with 
samarium  doped  ceria  (SDC)/scandia-stabilized  zirconia  (ScSZ) 
bilayer  as  electrolytes,  fabricated  by  PLD.  However,  in-depth 
analysis  of  the  polarization  characteristics  was  not  pursued  in 
that  work.  In  the  continuing  effort  to  optimize  the  fabrication 
process  and  to  improve  the  cell  performance,  a  comprehensive 
analysis  based  on  EIS  data  and  polarization  curves  is  given  in 
this  paper  over  the  temperature  range  of  400-600  °C.  Two  rep¬ 
resentative  metal-supported  SOFC  single  cells,  previous  Cell  #1 
[8]  and  recently  Cell  #2  fabricated  by  suspension  plasma  spray, 
were  chosen  as  examples  for  such  a  purpose. 

2.  Experimental  procedures 

Metal- supported  SOFC  single  cells  were  fabricated  using 
different  deposition  techniques.  Commercially  available  porous 
SS430  and  porous  Hastelloy  X  alloys  were  adopted  as  substrates 
for  selected  Cell  #1  and  Cell  #2,  respectively.  Porosity  measure¬ 
ments  of  the  SS430  and  Hastelloy  X  were  performed  using  the 
Archimedes  method  and  their  values  were  22  and  28%,  respec¬ 
tively.  Cell  #1  consists  of  NiO-SDC  (50:50  wt.%)  as  the  anode, 
ScSz  and  SDC  bi-layer  as  the  electrolyte,  and  Smo.sSro.sCoOs-^ 
(SSCo)-SDC  (75:25  wt.%)  as  the  cathode.  Cell  #2  consists  of 
NiO-SDC  (70:30  wt.%)  as  the  anode,  SDC  as  the  electrolyte, 
and  Smo.5Sro.5Co03_5  (SSCo)-SDC  (75:25  wt.%)  as  the  cath¬ 
ode.  The  active  area  is  0.23  cm2  for  Cell  #1  and  0.34  cm2  for 
Cell  #2,  respectively.  The  details  about  material,  structure,  and 
processing  methods  for  both  Cell  #1  and  Cell  #2  are  listed  in 
Table  1. 

Both  cells  were  mounted  on  an  alumina  tube  using  Ceram- 
abond  552  (Aremco).  Pt  meshes  were  used  as  both  anode  and 
cathode  current  collectors.  Both  cells  were  heated  to  650  °C  at 
2  °C  min- 1  and  dwelled  at  650  °C  for  5  h  for  the  purpose  of  fully 
reducing  the  anodes  by  gradually  raising  the  hydrogen  concen¬ 
tration  (maintaining  3%  water  stream).  Subsequently,  these  two 


Table  1 

Summary  of  materials,  structures,  and  processing  methods  for  Cell  #1  and  Cell 
#2 


Components 

Substrate 

Anode 

Electrolyte 

Cathode 

Cell  #1 

Material 

Processing 

Thickness 

SS430 

NiO-SDC 
Screen  printing 
90  pm 

ScSz/SDC 

PLDa 

2/20  pm 

SSCo-SDC 
Screen  printing 
25  pm 

Cell  #2 

Material 

Processing 

Thickness 

Hastelloy  X 

NiO-SDC 

SPSb 

30  pm 

SDC 

SPS 

30  pm 

SSCo-SDC 
Screen  printing 
45  pm 

a  Pulsed  laser  deposition. 
b  Suspension  plasma  spray. 


cells  were  heated  to  800  °C  at  3  °C  min-1  to  sinter  the  cathode 
for  30  min,  and  then  cooled  down  to  400  °C  at  2  °C  min-1 .  EIS 
and  current-voltage-power  curve  were  measured  from  600  to 
400  °C  at  interval  of  50  °C. 

Voltage-current  (/-V)  curves  were  measured  using  Solartron 
1480A  potentiostat  with  a  scan  rate  of  4  mV  s-1  in  the  poten¬ 
tial  range  from  open  circuit  voltage  (OCV)  to  0.3  V.  EIS  was 
measured  over  the  frequency  range  of  100  kHz  to  0.1  Hz  with  a 
Solartron  1260  frequency  response  analyzer  (FRA)  connected 
to  a  Solartron  1480 A  potentiostat.  The  perturbation  amplitude 
was  set  to  50  mV  at  open  circuit.  The  requirement  of  pertur¬ 
bation  amplitude  is  determined  by  the  internal  characteristics  of 
the  investigated  system,  to  a  great  extend,  by  the  linearity  degree 
of  the  investigated  system.  In  theory,  with  smaller  amplitude  the 
system  is  perturbed,  in  better  precision  the  linearization  process 
could  be  implemented.  In  engineering,  larger  amplitude  signal 
is  expected  because  perturbation  signal  with  smaller  amplitude 
corresponds  to  stricter  test  requirement  of  high-signal  noise 
ratio.  SOFCs  usually  exhibit  strong  linearity  over  a  large  current 
polarization  range.  Consequently,  signals  with  larger  amplitude 
than  the  traditional  10  mV  could  be  applied  to  impedance  mea¬ 
surement  for  SOFCs  with  tolerable  error.  For  each  designed 
temperature  of  600,  550,  500,  450  and  400  °C,  both  impedance 
and  I-V  data  were  measured  twice  at  10-min  interval.  All  above 
measurements  were  under  conditions  of  moist  hydrogen  as  the 
fuel  and  dry  air  as  the  oxidant. 

The  tested  cells  were  immersed  into  epoxy  followed  by  solid¬ 
ification,  and  then  sectioned  and  polished  with  different  diamond 
slurries.  The  cross-sectional  microstructures  of  two  cells  were 
then  examined  with  a  scanning  electron  microscope  (SEM, 
Hitachi  S-3500N)  at  two  different  resolutions  described  below. 

3.  Experimental  results  and  discussions 

3.1.  Cell  micro  structures 

In  order  to  investigate  microstructures  affected  by  PLD  and 
SPS,  the  two  cells  were  characterized  using  SEM  at  different 
resolutions. 

Fig.  1  clearly  shows  multi-layer  micro  structures  of  the  cross- 
sections  for  Cell  #1  and  Cell  #2  fabricated  in  lab  scale.  For  both 
cells,  the  electrode  layers  are  all  porous,  which  are  necessary  for 
reactant  gas  feeds  and  diffuses.  However,  the  interfaces  between 
layers  from  the  two  cells  are  drastically  different  due  to  the 
processing  methods. 

In  Fig.  la  and  b  for  Cell  #1,  the  AFL  is  much  finer,  not 
only  extending  the  length  of  triple  phase  boundary  to  speed 
up  the  absorption/adsorption  reaction  of  oxygen  ion,  but  also 
enhancing  the  bonding  of  electrolyte/anode  interface  to  reduce 
the  interfacial  contact  resistance.  Our  previous  work  indicated 
that  various  microstructures,  density  and  interfacial  properties 
of  bi-layer  ScSZ/SDC  film  could  be  achieved  by  controlling  the 
process  parameters  of  PLD  technique  [7].  The  bi-layer  elec¬ 
trolyte  ScSZ/SDC  fabricated  by  PLD  at  a  very  low  temperature 
is  fine  and  dense  without  any  pinholes  or  cracks,  which  not 
only  improves  OCV  [9],  but  also  meets  the  requirement  of  gas 
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Fig.  1.  Cross-sectional  SEM  images  of  the  two  tested  cells:  (a)  350x  magnification  for  Cell  #1;  (b)  2000x  magnification  for  Cell  #1;  (c)  150x  magnification  for 
Cell  #2;  (d)  lOOOx  magnification  for  Cell  #2. 


tightness.  From  the  SEM  images  in  Fig.  1 ,  it  can  be  seen  that  no 
visible  separation  gaps  between  interfaces  of  both  anode  side  and 
cathode  side  for  Cell  #1  can  be  found,  suggesting  the  interface 
contacts  are  excellent. 

Fig.  lc  and  d  shows  the  microstructure  of  the  cross-section 
for  Cell  #2  fabricated  by  SPS  method,  which  has  been  well 
proven  to  be  cost-effective  and  suitable  for  scale-up  and  offers  an 
alternative  technique  to  deposit  SOFC  functional  layers  on  metal 
substrate  [9].  Several  tiny  cracks  are  found  in  electrolyte  layers 
for  Cell  #2,  but  none  of  them  penetrated  completely  through  the 
overall  layer.  In  Fig.  Id,  no  clear  interfacial  gap  can  be  found  on 


both  sides  of  the  anode  and  the  cathode  for  Cell  #2,  suggesting 
relative  strong  interfacial  contact.  However,  compared  with  Cell 
#1,  the  interfacial  bonding  of  Cell  #2  seems  not  as  compact  as 
that  of  Cell  #1 .  These  interfacial  differences  could  result  in  large 
differences  in  ohmic  and  polarization  resistances  of  these  two 
cells.  This  will  be  further  discussed  in  Section  3.3. 

3.2.  Electrochemical  polarization 

Fig.  2  shows  the  voltage-current  and  power-current  curves 
for  both  Cell  #1  and  Cell  #2  over  the  temperature  range  of 
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Fig.  2.  Electrochemical  polarization  over  the  temperature  range  of  400-600  °C  for  Cell  #1  (a)  and  Cell  #2  (b). 
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Fig.  3.  (a)  Theoretical  and  measured  values  for  open  circuit  voltage,  where  OCV  is  for  the  theoretical  value,  OCV1  and  OCV2  for  the  measured  values  of  Cell  #1 
and  Cell  #2,  respectively  and  (b)  maximum  power  density  of  the  metal- supported  SOFCs  over  the  temperature  range  of  400-600  °C. 


400-600  °C.  A  thin  ScSZ  layer  was  employed  for  Cell  #1  to 
block  the  electronic  leaking  from  SDC  electrolyte  [8].  Cell  #1 
achieved  an  open  circuit  voltage  (OCV)  of 0.980-1 .024  V,  which 
is  much  higher  than  the  OCV  of  0.826-0.944  V  for  Cell  #2 
where  only  one  single  SDC  layer  of  electrolyte  was  used.  With 
respect  to  cell  performance  at  the  operation  temperatures  of  400, 
450,  500,  550  and  600  °C,  the  maximum  power  densities  of  Cell 
#1  are  0.063,  0.096,  0.124,  0.144,  and  0.160Wcm-2,  respec¬ 
tively.  The  maximum  power  densities  of  Cell  #2  are  0.016, 0.037, 
0.074,  0.123,  and  0.181  W  cm-2,  respectively. 

Fig.  3a  and  b  shows  the  theoretical  open  circuit  voltage  (OCV 
for  theoretical  value),  the  measured  OCV  (OCV  1  for  Cell  #1  and 
OCV2  for  Cell  #2),  and  the  measured  maximum  power  den¬ 
sity  (MPD)  related  to  operation  temperatures  for  the  two  cells, 
respectively.  As  expected,  the  measured  OCV  of  the  two  cells 
decrease  slightly  while  the  MPD  of  both  cells  increase  dramati¬ 
cally  with  the  increase  of  operation  temperature.  The  difference 
between  the  theoretical  values  and  the  experimental  ones  for 
the  OCV  could  be  contributed  from  different  sources,  such  as 
the  porosity  of  the  electrolyte  layer,  the  electronic  conductivity 


from  SDC  electrolyte,  or  the  fuel  leaking  from  the  edge  sealing. 
One  point  should  be  noted  that  MPD  of  Cell  #2  exceeds  that  of 
Cell  #1  at  600  °C  although  the  measured  OCV  of  Cell  #2  is  still 
lower  than  the  measured  OCV  of  Cell  #1 .  In  order  to  understand 
these  phenomena,  the  electrochemical  impedances  for  both  Cell 
#1  and  Cell  #2  were  measured  under  open  circuit  condition. 

3.3.  Electrochemical  impedance  measurement 

EIS  has  been  recognized  as  a  powerful  diagnostic  tool 
to  investigate  fundamental  processes  and  various  polarization 
losses  for  SOFCs  [20].  In  order  to  understand  and  compare  vari¬ 
ous  polarization  losses  and  associated  reaction  kinetics  for  both 
Cell  #1  and  Cell  #2,  electrochemical  impedances  were  mea¬ 
sured  under  open  circuit  condition  over  the  temperature  range 
of  400-600  °C.  Nyquist  plots  for  Cell  #1  and  Cell  #2  over  the 
temperature  range  of  400-600  °C  are  shown  in  Fig.  4a  and  b, 
respectively. 

Ohmic  resistances  R&  and  polarization  resistances  Rv  of 
the  two  cells  over  the  temperature  range  of  400-600  °C  were 


Fig.  4.  The  ac  impedance  spectra  of  two  metal- supported  cells  under  open  circuit  condition  over  the  temperature  range  of  400-600  °C  for  Cell  #1  (a)  and  Cell  #2  (b). 
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Temperature  (°C) 

Fig.  5.  (a)  Ohmic  resistances,  (b)  polarization  resistances  and  (c)  cell  resistances  of  the  two  metal- supported  cells  over  the  temperature  range  of  400-600  °C  under 
open  circuit  condition. 


extracted  from  the  measured  Nyquist  plots  in  Fig.  4,  assuming 
cell  resistance  Rt  =  Rp+R&.  Generally  speaking,  the  main  con¬ 
tributor  to  R&  is  from  the  electrolyte  [21] .  Fig.  5a-c  shows  curves 
of  R&,  Rp,  and  R{  as  a  function  of  the  operation  temperature  for 
the  two  cells,  respectively.  It  can  be  seen  that  both  of  RQ  and  Rp  of 
the  two  cells  decrease  with  the  increase  of  measurement  temper¬ 
ature.  Both  Rc  and  Rp  of  Cell  #2  are  obviously  greater  than  those 
of  Cell  #1  over  the  low  temperature  range  of  400-500  °C  likely 
due  to  Cell  #2  with  weaker  bonding  of  the  electrolyte-electrode 
interfaces  than  that  of  Cell  #1  as  shown  in  Fig.  1.  However,  RQ , 
Rp,  and  Rt  of  Cell  #2  are  more  sensitive  to  the  change  of  tem¬ 
perature  than  those  of  Cell  #1  due  to  the  energy  difference  of 
processing  between  PLD  and  SPS.  The  impact  from  fabrication 
processes  becomes  less  important  when  the  reaction  kinetics  get 
enhanced  at  elevated  temperature.  Consequently,  over  550  °C, 
the  large  difference  in  polarization  resistance  for  the  two  cells 
gradually  disappears  and  almost  converges  to  a  small  value.  In 
this  situation,  ohmic  loss  dominates  the  overpotential  losses.  At 
600  °C,  the  ohmic  resistance  of  Cell  #2  becomes  lower  than  that 
of  Cell  #1  and  this  is  most  likely  due  to  the  low-ionic  conductiv¬ 
ity  of  ScSZ  in  the  bi-layer  electrolytes  and  the  potential  reaction 
between  ScSz  and  SDC  in  Cell  #1  [7].  The  lower  cell  resistance 
of  Cell  #2  is  in  agreement  with  the  phenomenon  in  Section  3.2, 
where  MPD  of  Cell  #2  exceeds  that  of  Cell  #1  at  600  °C  although 
the  OCV  of  Cell  #2  is  still  lower  than  that  of  Cell  #1. 


3.4.  Modeling  of  cell  polarization 

3.4.1.  Analysis  of  polarization  characteristics 

In  the  simplest  case,  the  electrochemical  processes  charac¬ 
terized  for  Cell  #1  and  Cell  #2  in  Fig.  4  can  be  simulated 
by  an  equivalent  circuit  model  (ECM)  in  Fig.  6a  [22],  where 
(7?iCi)  simulating  high-frequency  processes,  ( R2C2 )  medium 
frequency  processes,  R0hm  ohmic  resistance  mainly  consisting 
of  electrolyte  resistance  and  various  contact  resistances,  L  an 


Fig.  6.  Two  equivalent  circuit  models  (ECM)  for  metal-supported  SOFC:  (a) 
ECM  to  simulate  the  electrochemical  processes  under  open  circuit  condition 
and  (b)  ECM  to  analyze  I-V  polarization  characteristics. 
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inductance  introduced  by  connecting  conductor.  The  time  con¬ 
stants  of  the  electrochemical  processes  are  given  by  RC  and 
can  be  determined  through  frequency  response  measurements 
[23].  It  is  also  well  known  that  I-V  characteristics  reflect  the 
overlapped  effect  of  ohmic  polarization,  concentration  polariza¬ 
tion,  and  activation  polarization.  Hence,  dc  measurement  and 
analysis  of  steady-state  I-V  characteristics  are  effective  ways 
to  investigate  kinetics  for  SOFCs  [21],  although  the  capaci¬ 
tive  and  inductance  parts  in  Fig.  6a  could  not  be  reflected  in 
I-V  characteristics.  The  ECM  in  Fig.  6b  is  applied  to  analyze 
I-V  polarization  characteristics  for  Cell  #1  and  Cell  #2  [24], 
where  E0 cv  open  circuit  voltage,  ^act  activation  polarization,  r]con 
concentration  polarization,  r]ohm  ohmic  polarization,  R\oad  load 
resistance,  V  cell  voltage,  i  current  density. 

Previous  experimental  results  [25,26]  suggested  the  overpo¬ 
tential  on  the  anode  side  was  extremely  small  in  comparison 
to  the  one  on  the  cathode  side.  Thus,  activation  polarization 
and  concentration  polarization  from  the  anode  side  could  be 
ignored  for  simplification  of  the  analysis.  As  a  result,  cell  voltage 
V(i)  —  Eocv  9 act  -  9 ohm  ~  9con,c  can  be  simplified. 


V  (i)  —  Eocv  ?7act,c  ^ohm  9con,c 


(1) 


In  this  work,  flow  rates  of  fuel  and  oxidant  were  held  high 
enough  compared  to  the  stoichiometry  required  by  the  corre¬ 
sponding  current  density,  thus  concentration  polarization  could 
be  simply  treated  either  as  a  constant  or  be  ignored.  This  allows 
further  cell  voltage  simplification: 


V  (i)  —  E0  cv  7?act,c  90hm 


(2) 


where  9'ohm  =  ^0hm  +  9con,c •  The  activation  polarization  for  the 
cathode  is  normally  expressed  by  the  well-known  Butler- Volmer 
equation  [23] 


i  =  io,c 


fin  Q,Fr]  c 
RT 


—  exp 


^  (1  /0^eT?7act,c  ^  ^ 

(3) 


where  i  is  current  density  (mAcm-2),  z'o,c  cathode  exchange 
current  density  (mAcm-2),  fi  the  transfer  coefficient  (usually 
0.5)  [27],  electrons  transferred  per  reaction  (here  2),  F  Faraday 
constant  (96,485  C  mol-1),  ^act,c  cathode  activation  polarization 
(V),  R  universal  gas  constant  (8.3 14  J  mol-1  k-1),  T temperature 
(K).  If  we  defined  90 (T)  as 

9o(T)  =  (4) 

finQF 


Then,  Eq.  (3)  can  be  expressed  as 


For  high-activation  polarization,  it  is  possible  that 


9act ,c  >  90 (T)  (6) 

In  such  case,  Eq.  (3)  can  be  simplified  as  the  logarithmic 
current-potential  relation. 


RT  RT 

^act,c  =  -- — -  In  10, c  +  - — -  In  1 
finQF  finQF 


(7) 


Eq.  (7)  is  the  well-known  Tafel  equation  for  the  cathode.  One 
point  has  to  be  noted  that  Eq.  (7)  is  only  valid  for  i  >  z'o,c  [28]. 
For  low-activation  polarization,  it  is  possible  that 


7? act,  c  <  *70  (T) 


(8) 


In  such  case,  Butler-Volmer  equation  can  be  simplified  as 
the  linear  current-potential  relation  in  terms  of  Taylor  series 
expansion. 

RT 

^7act  ,c  —  i  —  (9) 

nQFio,c 


By  differentiating  Eq.  (9),  we  can  define 


DOCV  _  def^act ,c 
KCt,C  -  n. 


RT 

n&Eio,c 


(10) 


where  R°f^  is  cathode  charge  transfer  resistance  under  open 
circuit  condition  [29,30].  When  both  concentration  polarization 
of  the  whole  cell  and  activation  polarization  on  the  anode  side 
are  neglected  as  discussed  above,  polarization  resistances  Rp 
obtained  from  EIS  measurement  under  open  circuit  condition 
for  Cell  #1  and  Cell  #2  can  be  approximately  regarded  as 
i.e. 


R° 


R n 


(11) 


Based  on  Eq.  (10)  and  the  approximation  of  via  Eq. 
(11),  the  cathode  exchange  current  density  z'o,c  could  be  cal¬ 
culated.  The  calculated  exchange  current  density  z'o,c  and  the 
current  density  i*  corresponding  to  MPD  for  Cell  #1  and  Cell  #2 
over  the  temperature  range  of  400-600  °C  are  listed  in  Table  2. 
It  is  obvious  that  cathode  exchange  current  density  z'o,c  increases 
dramatically  with  the  increase  in  operating  temperatures  for  both 
cells,  so  does  the  current  density  f  ,  and  the  z'o,c  for  Cell  #1  is 
greater  than  the  z'o,c  for  Cell  #2  at  given  temperatures.  Subse¬ 
quently,  the  increase  of  the  cathode  exchange  current  density  z'o,c 
as  a  function  of  temperature  extend  the  range  of  current  density 
that  is  not  valid  for  the  Tafel  equation,  but  is  still  valid  for  linear 
activation  polarization  via  Eq.  (9)  as  discussed  above.  Moreover, 
the  induced  cathode  exchange  current  density,  as  functions  of 
both  operation  temperatures  and  fabrication  processes,  offers  a 
viewpoint  to  reveal  the  polarization  related  to  operating  tem¬ 
peratures,  fabrication  processes,  and  materials  properties  for 
SOFCs. 


Table  2 

Calculated  exchange  current  density  for  cathode  and  measured  current  density 
corresponding  to  MPD  at  operation  temperatures  from  400  to  600  °C 

Temperature  (°C) 


400 

450 

500 

550 

600 

Cell  #1 

z'o,c  (mA  cm  2)a 

21.3 

53.0 

114.1 

190.2 

298.1 

1  (mAcm-2)b 

108.0 

173.0 

225.0 

263.0 

308.0 

Cell  #2 

io,c  (mA  cm  2) 

4.0 

12.4 

37.0 

92.5 

184.7 

f  (mAcm-2) 

26.0 

71.0 

153.0 

257.0 

414.0 

a  Cathode  exchange  current  density  (mA  cm  2). 
b  Current  density  corresponding  to  MPD  (mAcm-2). 
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Table  3 

Calculated  values  of  rjo(T),  ^*act  and  rfe act  at  °Peration  temperatures  of 
400-600 °C 


Temperature  (°C) 

400 

450 

500 

550 

600 

m(T)  (mV)a 

58.0 

62.3 

66.6 

70.9 

75.2 

^ac,,c(mV)b 

136.5 

99.2 

66.7 

48.3 

40.5 

173.3 

160.1 

133.6 

108.1 

102.1 

b  Cathode  activation  polarization  at  current  density  i  for  Cell  #1. 
c  Cathode  activation  polarization  at  current  density  i  for  Cell  #2. 


With  Eq.  (2),  we  can  get  rjacUc(i)  =  Eocv  -  V(i)  -  nohm(i), 
where  i  is  current  density  at  different  temperatures  of 
400-600  °C  as  listed  in  Table  2,  rfohm(i)  ~  i  Rohm,  R0hm  ohmic 
resistance  obtained  from  EIS  measurement.  Consequently,  we 
could  approximately  calculate  ^act,c(0  at  different  temperatures 
of  400-600  °C  in  order  to  investigate  the  cathode  activation 
polarization  over  the  current  density  range  of  0  ~  i  .  The  val¬ 
ues  of  X] o(7),  r]iact  and  r]2act  at  current  density  of  f  over  the 
temperature  range  of  400-600  °C  were  listed  in  Table  3.  The 
calculated  results  listed  in  Table  3  indicated  that  Eq.  (8)  was 
not  strictly  applicable  at  current  density  of  f .  Two  special  situ¬ 
ations  including  high- activation  polarization  and  low- activation 
polarization  have  been  discussed  previously  [27,31].  In  terms  of 
the  applicable  range  of  these  two  simplified  activation  models 
via  Eqs.  (7)  and  (9)  discussed  by  Chan  et  al.  [27],  the  low- 
activation  polarization  situation  is  more  suitable  to  the  calculated 
cathode  activation  polarization  r]*a  c  at  current  density  f .  On 
one  hand,  activation  polarization  corresponding  to  the  current 
density  f  is  maximized  in  the  range  of  0  ~  f .  On  the  other 
hand,  considering  the  activation  polarizations  on  the  anode  side 
and  concentration  polarization  of  the  cell  are  ignored,  the  cal¬ 
culated  cathode  activation  polarization  rjsict,c  is  expected  to  be 
greater  than  the  actual  activation  polarization.  As  a  result,  with 
respect  to  the  investigated  current  density  range  of  0~z*,  it 
is  reasonable  to  employ  linear  model  via  Eq.  (9)  to  approxi¬ 
mate  the  activation  polarization  of  Cell  #1  and  Cell  #2.  The 
reduction  of  the  activation  polarization  with  increase  of  tem¬ 
perature  strengthens  the  linearity  of  activation  polarization  as 
indicated  in  Eq.  (8).  Table  3  suggests  that  Cell  #1  exhib¬ 
ited  less  activation  polarization  than  Cell  #2  at  each  a  given 
temperature. 

3.4.2.  Fitting  of  polarization  curves 

According  to  the  discussion  for  the  activation  polarization 
above,  Eq.  (2)  can  be  completely  linearized  over  the  investigated 


current  range  of  0  ~  i  (mA). 

V(i)  —  £ocv  —  iRact,c  —  IR  ohm  =  ^ocv  —  ^t 

(12) 

where 

^t  =  Ract,c  4“  ^ohm 

(13) 

With  the  measured  current-voltage  data  and  approximated 
Eq.  (12),  the  fitted  values  of  parameters  such  as  open  circuit 


Table  4 

Fitted  E'ocw  and  measured  E0 cv  at  operation  temperatures  of  400-600  °C 


Temperature  (°C) 


400 

450 

500 

550 

600 

Cell  #1 

£ocva  (V) 

1.0244 

1.0238 

1.0167 

1.0014 

0.9801 

Kc7  (V) 

1.0241 

1.0247 

1.0180 

1.0023 

0.9813 

Cell  #2 

Eocv  (V) 

0.9440 

0.9291 

0.9040 

0.8701 

0.8258 

K cv  (V) 

0.9404 

0.9257 

0.9054 

0.8740 

0.8316 

a  Measured  OCV. 
b  Fitted  OCV  from  I-V  curves. 


voltage  E'ocy  and  cell  resistance  F't  could  be  obtained  through 
the  fitting  function  Polyfit  in  Matlab,  as  listed  in  Tables  4  and  5, 
respectively.  Subsequently,  the  normalized  and  accumulated- 
residuals  (NA-Res)  of  I-V  curve  fitting  with  Eq.  (12)  over  the 
temperature  range  of  400-600  °C  are  calculated  via 

100 

NA  —  Res  =  /  J  ^measured  —  ^calculated  It  (14) 

i 

where  ^measured  is  the  measured  value  of  cell  voltage,  VCaicuiated 
the  calculated  value  of  cell  voltage  via  Eq.  (12)  with  parameters 
from  I-V  curve  fitting,  N  the  number  of  data  points  for  current 
density  range  of  0  ~  f  and  T  is  the  temperature. 

Tables  4  and  5  indicated  that  the  fitted  cell  resistance  and 
open  circuit  voltage  agreed  well  with  the  measured  cell  resis¬ 
tance  and  open  circuit  voltage.  Table  6  shows  that  normalized 


Table  5 

Fitted  R[  and  measured  Rt  at  operation  temperatures  of  400-600  °C 


Temperature  (°C) 


400 

450 

500 

550 

600 

Cell  #1 

Rta  (Q  cm2) 

4.2854 

2.7360 

2.0784 

1.7459 

1.4951 

R'th  (Q  cm2) 

4.2273 

2.7348 

2.0856 

1.7498 

1.5025 

Cell  #2 

Rt  (£2  cm2) 

15.3496 

6.0884 

2.7851 

1.5042 

0.8961 

(£2  cm2) 

14.8022 

5.8597 

2.7918 

1.5564 

0.9511 

a  Cell  resistance  from  EIS  measurement  under  open  circuit  condition. 
b  Cell  resistance  from  I-V  fitting  via  Eq.  (12). 


Table  6 

Normalized  and  accumulated-residuals  (NA-Res)  of  I-V  fitting  at  operation 
temperatures  of  400-600  °C 

Temperature  (°C) 


400 

450 

500 

550 

600 

NA-Resa  (V) 
Cell  #1 

0.1334 

0.0735 

0.0555 

0.0462 

0.0398 

Cell  #2 

0.4052 

0.2195 

0.0932 

0.0855 

0.1406 

a  NA  -  Res 

II 

M 

1  ^measured  — 

^calculated  1 T  • 
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Table  7 

Measured  MPD  and  simulated  MPD'  at  different  operational  temperatures  of 
400-600 °C 


Temperature  (°C) 


400 

450 

500 

550 

600 

Cell  #1 

MPDa  (WcnU2) 

0.0626 

0.0963 

0.1244 

0.1435 

0.1600 

MPD'  b  (Won"2) 

0.0606 

0.0937 

0.1220 

0.1432 

0.1633 

Cell  #2 

MPD  (W  cm-2) 

0.0155 

0.0374 

0.0738 

0.1227 

0.1810 

MPD'  (W  cm-2) 

0.0159 

0.0395 

0.0811 

0.1404 

0.2186 

a  Measured  MPD. 
b  Simulated  MPD. 


and  accumulated-residuals  of  linear  fitting  are  quite  small  with 
ranges  from  0.0398  to  0.4052  V  related  to  each  100  current 
points  at  the  temperatures  of  400-600  °C.  The  fitting  residuals 
for  both  cells  decrease  with  the  increase  of  operating  tem¬ 
peratures.  This  trend  is  consistent  with  the  result  in  Table  2. 
The  only  exception  occurs  at  the  temperature  of  600  °C  for 
Cell  #2,  which  is  most  likely  due  to  the  instability  of  newly 
fabricated  metal-supported  cell.  Moreover,  the  normalized  and 
accumulated-residuals  NA-Res  of  Cell  #1  is  much  lower  than 
that  of  Cell  #2  at  each  given  temperature.  All  these  results 
strongly  suggest  that  the  approximation  and  discussion  on  the 
activation  polarization  and  I-V  characteristics  above  is  reason¬ 
able  and  feasible. 

With  the  fitted  open  circuit  voltage  E'ocy  and  cell  resistance 
through  fitting  at  temperatures  of  400-600  °C,  the  cell  MPD  can 
be  predicted  in  terms  of  Pmax  =  (l/2)(£QCV2/R()conveniently. 
As  indicated  in  Table  7,  Cell  #1  exhibits  higher  precision  in 
MPD  prediction  than  Cell  #2  at  given  operation  temperatures. 
This  result  is  also  consistent  with  the  conclusion  drawn  from 
Table  6  for  Cell  #1  and  Cell  #2  over  the  temperature  range  of 
400-600  °C.  The  above  analysis  of  polarization  characteristics 
related  to  operating  temperatures  and  fabrication  processes  is 
based  on  fitting  residuals  rather  than  the  intermediate  variables 
including  fitted  resistance  and  voltage,  or  by-production  such  as 
the  predicted  MPD. 

4.  Conclusions 

Polarization  characteristics  related  to  fabrication  processes 
and  operating  temperatures  have  been  investigated  for  metal- 
supported  SOFCs  fabricated  by  PLD  and  SPS,  respectively.  With 
EIS  and  current-voltage  polarization  curves,  the  electrochemi¬ 
cal  mechanism  for  the  linear  polarization  characteristics  of  the 
metal- supported  SOFCs  has  been  addressed  in  terms  of  cath¬ 
ode  exchange  current  density,  polarization  loss  and  maximum 
power  density  over  the  temperature  range  of 400-600  °C.  Exper¬ 
imental  and  simulating  results  indicated  that  high-operating 
temperatures  could  extend  the  range  of  current  density  with  lin¬ 
ear  activation  polarization,  and  low- activation  polarization  due 
to  fabrication  processes  could  also  strengthen  the  linearity  of 
activation  polarization.  In  summary,  fabrication  processes  and 
operating  temperatures  play  an  important  role  in  the  linearity 


of  the  polarization  characteristics  of  metal- supported  SOFCs. 
The  presented  method  in  this  manuscript  also  offers  a  viewpoint 
to  reveal  the  polarization  characteristics  of  SOFCs  related  to 
materials  properties,  fabrication  processes,  and  operation  tem¬ 
peratures. 
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